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ABSTRACT 

We examine the ability of the future Planck mission to provide a catalogue of galaxy 
clusters observed via their Sunyaev-Zel'dovich distortion in the cosmic microwave 
background. For this purpose we produce full-sky Sunyaev-Zel'dovich maps based on 
N-body simulations and scaling relations between cluster properties for several cos- 
mological models. We extrapolate the N-body simulations by a mass function to high 
redshifts in order to obtain a realistic SZ background. The simulated Planck observa- 
tions include, besides the thermal and kinematic Sunyaev-Zel'dovich effects, contribu- 
tions from the primordial CMB, extragalactic point sources as well as Galactic dust, 
free-free and synchrotron emission. A harmonic-space maximum-entropy method is 
used to separate the SZ signal from contaminating components in combination with a 
cluster detection algorithm based on thresholding and flux integration to identify clus- 
ters and to obtain their fluxes. We estimate a survey sensitivity limit (depending on 
the quality of the recovered cluster flux) and provide cluster survey completeness and 
purity estimates. We find that given our modelling and detection algorithm Planck 
will reliably detect at least several thousands of clusters over the full sky. The exact 
number depends on the particular cosmological model (up to 10000 cluster detections 
in a concordance ACDM model with as = 0.9). We show that the Galaxy does not 
significantly affect the cluster detection. Furthermore, the dependence of the thermal 
SZ power spectrum on the matter variance on scales of 8 h~^Mpc and the quality of its 
reconstruction by the employed method are investigated. Our simulations suggest that 
the Planck cluster sample will not only be useful as a basis for follow-up observations, 
but also will have the ability to provide constraints on cosmological parameters. 

Key words: cosmology: large-scale structure of the Universe - cosmology: cosmic 
microwave background - cosmology: theory - methods: data analysis - methods: sta- 
tistical. 



1 INTRODUCTION 

In addition to measuring primordial CMB fluctuations, the 
Planck Surveyor satellite, which is scheduled for launch 
in 2007, will have the ability to construct a cluster cat- 
alogue via the Sunyaev-Zel'dovich (SZ) effect (Sunyaev & 
Zel'dovich 1972; Sunyaev & Zel'dovich 1980; recent reviews: 
Rephaeli 1995; Birkinshaw 1999; Carlstrom, Holder, & Reese 
2002). Since Planck will provide detailed full-sky maps at 
nine different frequencies, ranging from 30 to 850 GHz, it is 
a particularly suitable instrument to detect the thermal SZ 
effect of galaxy clusters owing to its distinctive frequency 
dependence. Moreover, since the SZ effect is redshift inde- 
pendent and clusters - as the most massive collapsed objects 
in the universe - are tracers of the underlying linear matter 
density perturbations, their observation provides a unique 
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tool for studying large scale structure and constraining cos- 
mological parameters. 

In this work we investigate the yield, completeness and 
purity of such a Planck SZ survey for a range of cosmologi- 
cal models and estimate the cosmological significance of such 
the resulting cluster catalogue. Furthermore, based on our 
simulations and cluster selection we estimate a Planck clus- 
ter survey flux limit and provide tabulated cluster number 
counts for each model. In contrast to previous work (Kay, 
Liddle, k Thomas 2001; Diego et al. 2002; White 2003), 
which studied the cluster sample either on sky patches or 
neglected contaminants, we perform full-sky simulations in- 
cluding all known significant extragalactic and Galactic con- 
taminating components. Our modelling of the SZ effect as 
well as contaminants aims to be as realistic as current ob- 
servational constraints allow. To recover the SZ signal of 
single clusters we make use of the harmonic-space maxi- 
mum entropy method introduced by Stolyarov et al. (2002) 
in combination with a cluster finding algorithm based on 
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peak detection, thresholding and flux integration. As Stol- 
yarov et al. (2002) show the HSMEM is a novel method to 
separate various CMB components. The cluster finding al- 
gorithm is shown in this paper to be able to obtain reliable 
flux as well as radius estimates for a cluster sample above the 
estimated survey flux limit. Note that besides the assumed 
cosmological model the cluster recovery algorithm has ma- 
jor implications on the number of clusters obtained by the 
Planck survey and thus its impact on constraining cosmolog- 
ical parameters. It is therefore important to explore, improve 
and combine techniques in order to find the optimal method 
and maximise the survey yield. Since our modelling exhibits 
a high degree of realism, the given estimate of the cluster 
number which can be obtained by the Planck survey using 
our algorithm is reliable. Even though further improvement 
is possible, the presented combined cluster detection method 
- as we show below - provides already the quality to draw 
interesting cosmological conclusions from the obtained clus- 
ter sample. 

The outline of the paper is as follows. In section |21 a 
short introduction to the SZ effect is given. The modelling 
of full-sky SZ maps of different cosmological models and of 
the contaminants is described in sectional Moreover, section 
Elcontains a brief description of how Planck observation sim- 
ulations are obtained from the modelled components. In sec- 
tion^l the applied reconstruction method used to recover the 
fluxes of single galaxy clusters is explained. Section |5] con- 
tains a presentation and discussion of the results. Finally, a 
summary and conclusions are given in section |S] 



Model ^0 ^h 



as m/10i2/i-iM(; 



rCDM 


1.0 


0.0 


0.5 


0.21 


0.6 


2.22 


ACDM 


0.3 


0.7 


0.7 


0.21 


0.9 


2.25 



Table 1. Key parameters of the Hubble Volume simulations. r2o 
is the density parameter, Qa = A/3Hq the cosmological con- 
stant, h = ifo/lOOkms"-*^ Mpc"-*^ the Hubble constant, F the 
power spectrum shape parameter, ag the variance of linear fluc- 
tuations smoothed with a top-hat filter on the (comoving) scale 
of 8 h~^ Mpc at 2; = and m the N-body particle mass. 



electrons has been estimated. The spectral shape of the sec- 
ond contribution in equation (0, the kinematic SZ effect, is 
given by 



h{x) 



(e^-l)2' 



(4) 



and its magnitude, f3 = t'pec/c, depends on the uniform 
peculiar line-of-sight bulk motion of the cluster's electron 
plasma, '^pec- 



: (JT 



Tie dl , 



(5) 



is the Thomson optical depth. In the case the cluster can be 
assumed to be isothermal the Comptonization parameter 
can be expressed by 



y 



KrrieC^ 



(6) 



In this paper we concentrate on the thermal SZ effect and 
treat the kinematic merely as a contaminant to the thermal 
SZ. 



2 THE SUNYAEV-ZEL'DOVICH EFFECT 

The CMB anisotropy caused by the SZ effect can be sepa- 
rated into two contributions which are distinguished by the 
origin of energy of the scattering electrons that is responsible 
for the shift of photon frequency. The total CMB distortion 
due to the SZ effect is given by 



—r- = g{x)y-h{x)(3r, 



(1) 



where x = hv/k^TQ with To = 2.725 K (Mather et al. 1999) 
and /o = 2k%T§ /h^c? . The first term in equation (^ is the 
so called thermal SZ effect due to the thermal motion of 
electrons of the intra- cluster gas. The thermal SZ effect has 
a spectral shape given by 



9{x) = 



(e^ 



e" + l 



(2) 



and a frequency independent magnitude, the Comptoniza- 
tion parameter. 



y- 



kBCTT 



I 



UeTe dl . 



(3) 



In hot clusters (Te > 5keV) the relativistic electrons present 
slightly modify the spectral shape of the thermal SZ effect 
(Challinor & Fasenby 1998). This resulting relativistic cor- 
rection has not been taken into account in this work, since its 
effect on the results presented is negligible. The detect abil- 
ity of the effect from thermal (e.g. Pointecouteau et al. 1998) 
and non-thermal (e.g. Ensslin & Hansen 2004) relativistic 



3 SIMULATING THE SUNYAEV-ZEL'DOVICH 
EFFECT ON THE FULL SKY 

3.1 The cluster halo distribution 

There are two possible ways to obtain the distribution of 
cluster halos in redshift space. The first one is based on as- 
suming a cluster mass function (see e.g. Press & Schechter 
(1974); Sheth & Tormen (1999); Jenkins et al. (2001); 
Evrard et al. (2002) for cluster mass functions). This ap- 
proach determines the number of clusters within a given 
mass range at redshift z per comoving volume element. 



dN 
dz 



AQ (z) 

dzdQ^^ 






dn{M, z) 
dM 



dM, 



(7) 



where dV /{dzdO) is the comoving volume element and 
{dn/dM)dM is the comoving density of clusters of mass M. 
The second approach is based on performing N-body 
simulations and extracting halos from them to obtain halo 
catalogues which include halo masses, three-dimensional po- 
sitions and the halo three-dimensional velocities. Such cat- 
alogues have been obtained by applying a spherical over- 
density algorithm for the 'Hubble volume' simulations (see 
Frenk et al. 2000), the largest N-body simulations yet per- 
formed. Since, besides the halo redshift distribution, N-body 
simulations also provide information about halo clustering 
and dynamics, we decided to utilise these N-body simula- 
tions. These simulations investigate structure formation in 
two cosmologies, a rCDM and ACDM model. They provide 
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Figure 1. The cosmology dependent differential surface cluster 
number count per redshift bin Az = 1 for different lower cluster 
mass limits. The upper panel shows the number counts obtained 
from the fiducial ACDM simulation and the lower panel the ones 
for fiducial rCDM model (see TableQ. A typical cluster detected 
by Planck has a mass of approximately 5xl0^'^h~^ M©. The lower 
mass clusters provide a realistic cluster-cluster confusion noise. 



halo catalogues of light cone datasets for a sphere, two oc- 
tants and a wedge for each cosmological model (see Table 
Q). To obtain a full-sky realisation of the SZ effect we make 
use of the sphere and octant datasets. Nevertheless, espe- 
cially in the case of the ACDM model, the redshift cut-off 
of even the octants (z ~ 1.5) is rather low, since objects of 
intermediate masses (M ?^ 1 x 10^^h~^ M© and below) ex- 
ist at higher redshift. These objects contribute significantly 
to the background cluster confusion, given the ACDM clus- 
ter redshift distribution (see Figure 0. Therefore we extend 
the N-body simulations to higher redshifts by using a mass 
function to obtain the differential number density, 



dn 
d\nM 



{M,z) 



Pmd\na{M) 



M 
exp{- 



X 



c^lnM 

B-\n[g{z)aiM)]n 



(8) 



whose calibration has been directly derived from the Hub- 
ble volume (HV) simulations (Evrard et al. 2002) and is 
thus consistent. These extra high redshift halos due to 
the extension of the simulations are mainly added to pro- 
vide a realistic SZ background confusion. In equation (|H|, 
pm = Qm{0)3HQ / (SttG) is the present mass density, a{M) 
the variance of the density field smoothed with a spherical 
top-hat filter enclosing mass M at the mean density and 
g{z) is the linear growth factor normalised to ^(0) = 1. In 
the case of the ACDM model Evrard et al. (2002) find the 
best fit mass function parameters to the HV simulations to 
he A = 0.22, B = 0.73 and e = 3.86 and for the rCDM 
cosmology A = 0.27, B = 0.65 and e = 3.77. 



3.2 Halo properties from scaling relations 

Besides the mass of the dark matter halo, which is directly 
obtained from the N-body simulations, further cluster prop- 
erties, such as the electron gas temperature Te and virial 
radius r^, are needed in order to model the SZ effect of 
galaxy clusters. 



3.2.1 The mass-temperature relation 

Assuming that clusters are self-similar, virialised and 
isothermal, the mass-temperature (M-T) relation can be de- 
rived analytically (Kaiser 1986) and takes the form: 



Me] 



1015/i-iM 







P 



{l + z){QmAc/nm{z)y/^keY 



3/2 



(9) 



where Md is the mass of the cluster within its virial radius, 
Ac is the mean overdensity inside the virial radius in units 
of critical density and f3 parametrises the scaling relation. 
In general, p may be a function of Te as well as cosmology 
and redshift. In the case of a flat cosmology (Qa + Qm = 1) 



187r^ + 82{Qm{z) - 1) - 39{Qm{z) - if 



(10) 



(Eke et al. 1996). While in theory f3 ^ 0.75^ and is there- 
fore independent of Te and cosmology, this result, which has 
also been found by numerical simulations (see Pierpaoli et 
al. 2001, for a listing of results), differs from the ones ob- 
tained by observations which favour lower normalisations, 
i.e. (3 ^ 0.54(Tew/6keV)°-^^, where Tew is the emission- 
weighted X-ray temperature of the intracluster gas. This ob- 
servational normalisation has been derived from Finoguenov 
et al. (2001) by converting the mass at a mean overdensity of 
500 times the critical density to the virial cluster mass (see 
Hu & Kravtsov 2003). Note that the findings of Finoguenov 
et al. (2001) are in good agreement with Allen et al. (2001). 
To do the conversion, it has been assumed that the dark mat- 
ter distribution, which constitutes most of the cluster mass, 
obeys a NEW profile, since CDM is the major cluster mass 
component and simulations show that its distribution is well 
described by the NEW profile (see e.g. Navarro et al. 1997). 
Although there has been some controversy, this assumption 
is supported by observations. Eor example, by testing differ- 
ent dark matter models against mass profiles derived from an 
XMM-Newton observation of the relaxed cluster A478 under 
the hypothesis of hydrostatic equilibrium, Pointecouteau et 
al. (2004) find that the NEW profile is significantly preferred 
to other models. To distribute intracluster electron gas, how- 
ever, we use a /3-model. This assumption is supported by 
X-ray observations of relaxed galaxy clusters (see e.g. Et- 
tori et al. 2004). In the case of isothermality and hydrostatic 
equilibrium it has been shown that the expected gas profile 
of a NEW dark matter halo is well approximated by such a 
/3-profile (Eke, Navarro & Erenk 1998). Note further, that in 
the case of isothermality, Te = Tew , the observational result 
predicts a lower mass at fixed Te than is predicted by the- 
ory. The assumption that clusters are isothermal is in good 
agreement with recent XMM-Newton observations of outer 



^ Note that the theoretical prediction is based on a number of 
assumptions, such as hydrostatic isothermality of the cluster and 
y^-profile adherence of the cluster gas (see e.g. Sarazin 1986). 
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cluster regions which find that the cluster temperature pro- 
files are isothermal within ±10 per cent up to approximately 
half the virial radius (Arnaud et al. 2003). 

In the following we will use both normalisations of the 
M-T relation, the observed as well as the theoretical one, to 
analyse the difference it makes and if a Planck SZ survey 
will be able to distinguish between them. 



3.2.2 The integrated Comptonization parameter Y and 
the total cluster flux density Su 

After specifying the cluster abundance and the mass- 
temperature relation, we relate the observable flux due to 
the thermal SZ effect to the cluster mass. Since the spectral 
behaviour of the thermal SZ effect is the same for all clusters 
(see equation |3>, the total SZ flux density can be expressed 
in terms of the integrated Comptonization parameter 



Y 



ydA, 



(11) 



where dQ = dA/r'^, rd is the angular diameter distance to 
the galaxy cluster and one integrates over the cluster's pro- 
jected area. Then the total flux density Su at frequency z^, 
which is usually given in units of Jy, is 



SJx) 



(e^ - 1)2 



e^ + 1 



Y = gix)Y, 



(12) 



where again x = hv/ksTQ. 

In our case, where we assume the cluster to be isother- 
mal, the integrated Comptonization parameter is further 
given by 



Y 



_ J2 '^d ^e 



:/(/ 



riedl I dA 



kBCFT 



'TeNe 



(13) 



where we have used equation (|H|, and A^e the total number 
of thermal electrons within the cluster. The total number of 
electrons within a cluster's virial radius is proportional to 
the virial mass of the cluster, 



Ne= { ^-^j^ ] fgasMcl, 



(14) 



where fgas is the baryonic gas mass fraction of the total 
cluster mass, fn is the hydrogen fraction of the baryonic 
mass - in the following we assume /h = 0.76 - and rup is 
the proton mass. Therefore, the integrated Comptonization 
parameter Y, or respectively the total flux density, of the 
cluster can be related to the cluster's virial mass by 



Y 



rrieC^ 



l + fH 

2mv 



fgas Mel 



0.1x(l + /.)/,..(^) 



Me. 



10i5/i-iM( 



X 



Td 



h arcmin 



100/i-iMpc^ 

0.1 x{l+fH) fgas (3-\l + z) 

1/3 / ^^ ^ 5/3 

X 



Qm{z) 



M, 



cl 



1015/i-iMc 



rd 



100/i-iMpc 



h arcmin 



(15) 



Hence, apart from the normalisation, Y depends on the clus- 
ter gas mass, the cosmology and the cluster redshift. 



3.2.3 The cluster electron gas density profile 

Since the /3-model has been empirically found to fit well the 
electron gas density profile obtained from X-ray observa- 
tions, and since most of the clusters will be unresolved given 
the resolutions of the Planck channels (see sect ion 15. 5. 6|) . we 
utilise a spherical isothermal /3-profile, 



Pir) = po 



1 + 



ii:)- 



-3/2/3p„f 



(16) 



to model the electron gas distribution inside the cluster's 
virial radius. The line-of-sight projected ^^-profile is given 
by: 



'(r) 



■ rieo 



1 + 



iB' 



l/2-3/2/3p„f 



(17) 



where re is the cluster core radius and Ueo is the central 
electron gas density. For our simulations we make the com- 
mon choice of /3prof = 2/3. While this choice of P is consis- 
tent with the assumption of hydrostatic equilibrium for the 
virial M-T relation normalisation, for the one derived from 
X-ray observations the resulting cluster ^-profile given the 
gas density mimics a temperature gradient. When a trunca- 
tion of the cluster electron gas profile (|16f) is performed at 
the virial radius, the projected profile becomes 



i(r) = 



2 He 



Vl + (r/rcf 



X arctani 



Ib-^-jr/rcf 
1 + (r/rcf 



(18) 



where b relates the cluster core radius to the virial radius 
(see below). This truncation ensures that the cluster SZ sig- 
nal drops to zero within a finite range from the cluster cen- 
tre. Thus, in the cluster outskirts the projected profile (fTS|l 
declines more rapidly than the one given by (fT7)) . The result- 
ing change in the shape of the projected profile also agrees 
better with predictions obtained from hydrodynamical sim- 
ulations. 

As mentioned above, we use the virial radius to define 
the edge of the cluster. Since Mei = fvrpr^, where p is the 
mean matter density inside r^ , the virial radius of the cluster 
is given by: 



2.12 X 



Mel 



1/3 



HmA, 



^m{z) 



W^H-^Mq 

1/3 

h~^Mpc. 



(1 + ^)"' 



(19) 



Moreover, we relate the cluster core radius re to its virial 
radius rv by re{z) = b{z) rv, where b{z = 0) ^ 0.14 and the 
evolution of the cluster profile can be described by: 



,(z)^re{z^O) X (1 + z)'/^ 



(20) 



The scaling b{z) has been chosen so that the evolution of 
the central electron density of the cluster Ueo is consistent 
with simulations (see e.g. Eke et al. (1998)). 
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Given the shape of the electron gas density profile, the 
electron gas distribution and therefore the optical depth can 
be normalised in terms of the cluster's virial mass by utilis- 
ing equations (|14|) and (|16|) and by making use of 



f 

Jo 



gas{r)47Tr^ dr = fgasMd. 



(21) 



Hence, one deduces that the central electron gas density is 
given by: 



neo = 1.6 X lO'/i'/pas (^^/^) ( -^ 



Mpc 



10i5/i-iM( 
where w is given by 
w — b~ + arctan [h) — 



TV 

2 



(22) 



(23) 



when we assume a /3-profile with y^prof = 2/3. Note that - 
besides the relation for the central electron gas density - all 
scaling relations between cluster properties, which are given 
above, are independent of the chosen electron gas profile. 

Since the scaling relations presented above are for virial 
cluster properties, we have to rescale the cluster masses 
given by the Hubble volume cluster catalogues and the mass 
function of Evrard et al. (2002), as they provide the cluster 
mass within a radius for which the mean density contrast 
Ac = 200. The cluster mass rescaling is discussed in section 

El 

3.3 Varying ag 

In order to rescale the cluster halo dark matter masses to 
values of as which differ from the default as used in the 
conducted HV N-body simulations we follow the formalism 
presented in Evrard et al. (2002). A relation between the 
change in as, Aas, and the related change in cluster mass, 
AM, can be derived by setting to zero the total derivative of 
the differential comoving cluster number density (equation 
IHI at fixed redshift z. The resulting relation is given by 



Aag ,.,.. AM 

= a (M) — -— 

(78 ^ ^ M 

with 



a{M) 



aeff(M) + - 



(l-2V«eff(M)) 



(24) 



(25) 



e(lncr-i(M) + B)^-i 

where aeff(M) = diner ~^(M)/(ilnM and can be obtained 
by employing the quadratic relation describing the filtered 
power spectrum shape given in Evrard et al. (2002). As has 
been pointed out by Evrard et al. (2002), the number of the 
rarest most massive clusters is most sensitive to a change in 
as. 

This expression (equation I25f) differs from the one in 
Evrard et al. (2002), since we do not make any approxima- 
tion. In particular, for lower masses this expression is slightly 
more accurate (10 per cent for Md = 1 x 10^^ /i~^ M©). The 
presented rescaling relation between the halo dark matter 
mass and the change in as has been used to convert the 
CDM cluster masses of the fiducial HV simulations (see Ta- 
ble Q) to the ones of models with different values of as , a 
still rather unconstrained parameter by today's observation. 



3.4 Converting to virial mass 

Commonly, a cluster halo is defined to be a spherical over- 
dense region - in the case of N-body simulations - of par- 
ticles, whose internal mean density, (p), is given as a mul- 
tiple either of the critical density pc{z) — 3H(z)^ /SttG or 
of the overall mean matter density p(z) = Qrn{z) pc{z). For 
the HV catalogues Ac = 200, where {p) = AcPc(^)- In the 
case of the rCDM model this is not far-off the prediction 
of Ac ~ 178 given by the spherical top-hat collapse model 
for a sphere that encloses the virial mass, Md (see Lacey 
& Cole (1993)). For the fiducial ACDM cosmology however 
this value differs approximately by a factor of two from the 
predicted collapse model value of Ac(0) ^ 103 at redshift 
2; = 0, and Ap(0) ^ 667 for the HV simulations. Since the 
scaling relations presented above, especially the M-T rela- 
tion (see equation 13, are derived for the virial cluster prop- 
erties, either the cluster mass has to be directly converted to 
yield the virial mass of the cluster or one has to rescale the 
relations by deriving an approximate formula to the cluster 
mass conversion. Note, in either case a halo density profile 
has to be assumed in order to be able to convert various 
definitions of the cluster mass into each other. 

Therefore, in this paper we proceed by assuming a NEW 
halo profile of the dark matter which contributes mainly 
to the total cluster mass and follow the recipe as it has 
been suggested in Hu & Kravtsov (2003) to rescale the dark 
matter halo masses. 



3.5 Simulating Planck Observations 

The thermal and kinematic SZ effect contributions to a full- 
sky map at each Planck channel are obtained from cata- 
logues of cluster properties by employing the relations and 
profile discussed above. Further components regarded as 
contamination of the clusters' SZ signal were also simulated. 
The following paragraphs contain short descriptions of the 
realisation of these further components. 



3.5.1 Primordial CMB 

The realizations of full-sky maps of the primordial CMB 
for different cosmological models are generated by assuming 
that the primordial CMB is a homogeneous random Gaus- 
sian field entirely described by a standard CDM power spec- 
trum whose shape and normalisation agrees with WMAP 
measurements (Bennett et al. 2003). The Ci coefficients have 
been created using CAMB (Lewis, Challinor & Lasenby 
2000). 



3.5.2 Point sources 

Point source contamination has been included in our sim- 
ulations. Here we follow the modelling approach presented 
in Hobson et al. (1999). In respect of the component recon- 
struction, which is performed in spherical harmonics space 
(see section 14. If) , we make the assumption that the point 
source contribution in each Planck frequency channel should 
be well modelled as random Gaussian noise in the Fourier or 
spherical harmonics domain. The contribution by the point 
source population is determined by number counts predicted 
by Toffolatti et al. (1998) based on the evolution model of 
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Instrument 




LFI 










HFI 






Central frequency (GHz) 


30 


44 


70 


100 


143 


217 


353 


545 


857 


Angular resolution 
(FWHM, arcmin) 


33 


24 


14 


9.2 


7.1 


5 


5 


5 


5 


Sensitivity per pixel {fiK) 
(after three full sky coverages) 


4.5 


6 


10.5 


4.5 


4.9 


10.7 


32.7 


327.7 


14935 



Table 2. Basic observational characteristics of the 9 frequency channels of the Planck instruments (Low Frequency Instrument (LFI) 
and High Frequency Instrument (HFI)). The angular resolutions are quoted as FWHM for a Gaussian beam. The sensitivities quoted 
will be achieved by Planck after three full sky coverages. A pixel is a square whose side is the FWHM extent of the beam. 



radio sources of Danese et al. (1987) and the far-IR source 
evolution model by Franceschini et al. (1994), which has 
been updated by Burigana et al. (1997) to account for an 
isotropic sub-mm component. The rms thermodynamic tem- 
perature fluctuation within a defined pixel area due to the 
point source population in each Planck channel after beam- 
convolution is shown in Figure 4 of Hobson et al. (1999) for 
sources in the 10~^Jy< S{i') < 10 Jy flux range. We assume 
no correlation of the point sources with each other or with 
clusters of galaxies. In the case of radio point sources how- 
ever it is known that they have in fact a correlation with 
galaxy clusters. Thus, our assumption leads to some under- 
estimate of confusion at the lowest Planck frequencies. If IR 
sources originate at high redshift {z > 3) as suggested by 
current models as well as SCUBA observations at 850-/im 
(Holland et al. 1998; Smafl et al. 1998; Plain et al. 1998), the 
assumption of IR sources being uncorrelated with clusters 
which are detectable by Planck (see section 15. 2f) is valid. 
Generally, the contribution of clustered sources is found to 
be small in comparison to the Poisson term. Note that the 
Toffolatti model has been found to agree within a factor of 
two with results obtained from the WMAP mission (Ben- 
net et al. 2003). The WMAP results suggest that the model 
actually overpredicts the point source contribution. This is 
also supported by recent VSA 9C survey counts (see Fig- 
ure 13 in Bennett et al. 2003). Thus, our modelling is rather 
pessimistic concerning the radio point source contribution. 
IR sources are expected to dominate the number counts of 
bright objects at frequencies zy > 100 GHz, whereas radio 
sources are a significant contaminant below zy = 100 GHz. 
Currently the best observational limits on the IR source 
contribution are obtained from SCUBA. These observations 
suggest that the IR source contribution could be up to a 
few times higher than predicted by the Toffollati model at 
353 GHz. However, since the SCUBA surveyed area is rather 
small (a few 10 arcmin^) and IR sources relevant to Planck 
are at the bright end of the SCUBA population, the sample 
variance is quite large and does not yield any final conclu- 
sions about the IR source contribution. 



3.5.3 Galactic dust emission 

The dust emission is modelled using the DIRBE-IRAS 100- 
/im dust map (Schlegel, Finkbeiner & Davis 1998). Its res- 
olution of 5 arcmin is sufficient for our modelling purpose 
since it is comparable to the high frequency channel reso- 
lutions of Planck. To extrapolate the fiux at 100-/im to the 
Planck observing frequencies a one-component dust model 



with a temperature of Tdust = 18K and a dust emissivity 
;5 = 2 is assumed by taking into account the colour correc- 
tion of the DIRBE 100-/im filter. 



3.5.4 Galactic synchrotron emission 

The destriped version of the 408MHz Haslam survey 
(Haslam et al. 1982), to which artificial substructure on sub- 
degree scale has been added by extrapolating the power 
spectrum, is used as the synchrotron emission template. Its 
extrapolation to 300 GHz has been obtained by using an 
all-sky spectral index map derived from low frequency sur- 
veys at 408 MHz, 1420 MHz (Reich & Reich 1986) and 2326 
MHz (Jonas, Baart & Nicholson 1998) and padding the un- 
observed area around the South pole with the mean spectral 
index. By assuming a constant spectral index /3 = 0.9 the 
synchrotron emission at the Planck channel frequencies has 
been predicted. 



3. 5. 5 Galactic free-free emission 

A full-sky free-free emission template has been created on 
the basis of the DIRBE/IRAS dust map. 60 per cent of the 
free- free emission is assumed to be dust correlated. The dust 
uncorrelated free-free component has been obtained by fiip- 
ping the dust map north-south in Galactic coordinates. A 
spectral index of /3 = —0.16 has been assumed and the nor- 
malisation has been taken from Bouchet & Gispert (1999). 



3.5.6 Planck channel observations 

Realistic Planck surveyor simulations are then obtained at 
each of the n/ = 9 observing frequencies (see Table |2J 
along the Une-of-sight x by summing up the contributions 
of the physical components along this direction, convolving 
the resulting anisotropy map with the Planck beam at each 
particular frequency and adding the channel noise, which 
is assumed to be uncorrelated and Gaussian. This can be 
conveniently represented as a data vector of length n/ con- 
taining the resulting observing signal in direction x at each 
frequency: 



x) = Bty{x • x) y^^FtypSp{x) dQ' + ety{x), 



(26) 



where Fi^p = J tu(iy')fp{iy') dv' is the frequency response 
matrix, B^ is the beam profile of the z^th frequency channel. 
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Up-to-date beam sizes and instrumental channel noises for 
Planck are given on the Planck homepage^. 

In order to reduce the instrumental noise contribution 
in the high-frequency Planck channels to a level comparable 
to the point source noise, the presented results have been 
obtained by simulating Planck observations of one and a 
half year duration (approximately three full sky coverages) . 

Even more realistically one expects anisotropic instru- 
mental noise and spatial variations of the Galactic dust tem- 
perature and spectral index. Stolyarov et al. (2004) have re- 
cently extended their modelling to accommodate these ef- 
fects and applied the HSMEM to the simulated Planck chan- 
nel observations. Since they find that these spatial varia- 
tions do not affect the recovery quality of the components, 
our modelling is sufficient to obtain reliable results. In any 
case, these variations happen on scales larger than the typi- 
cal cluster scale and thus the assumption of locally constant 
noise and foreground properties in cluster-sized regions is a 
reasonable one. 



4 RECONSTRUCTING CLUSTER'S SZ 
SIGNAL 

4.1 The HSMEM 

In order to extract the thermal SZ signal from other CMB 
components on the sky, we use the harmonic-space maxi- 
mum entropy method (HSMEM) as presented in Stolyarov 
et al. (2002). This method allows us to perform a full-sky 
separation of all input components described above in spher- 
ical harmonic space. While the HSMEM has been developed 
mainly focusing on the recovery of the primordial CMB com- 
ponent and its power spectrum, it has been applied in this 
work to extract the thermal SZ component and all other 
components have been considered as contaminating noise. 
In the following a short description of the method is pre- 
sented. A detailed explanation of the method can be found 
in Stolyarov et al. (2002) and Hobson et al. (1998). 

The HSMEM is a Bayesian method. In the spheri- 
cal harmonic representation when we neglect any coupling 
between different (/,m) modes, the posterior distribution, 
which we want to explore, is given by 



Pr(a^ 



oc Pr(D£rn|a£m) Pr(a£m), 



(27) 



where Pr(D£rn|a€m) is the likelihood of the signal vector 
a^rn given the data D^rn- The prior probability Pr(a£rn) con- 
tains our previous knowledge or assumptions about the sig- 
nal vector before acquiring any data. In the case of the 
maximum-entropy method the specific feature that distin- 
guishes it from other Bayesian component separation algo- 
rithms is that it assumes the prior to be of entropic form, 
Pr(h£rn) oc exp[a5'(h£rn, m)], where the 'hidden' vector h^rn 
is related to the signal vector by a^rn = ^i^im and m is a 
model vector to which him defaults in the absence of data. 
5'(h£rn, in), whose form has been derived in Gull & Skilling 
(1990) and Hobson & Lasenby (1998), is the cross-entropy 
of him and m and a is a regularisation parameter, which is 
treated as another optimisation hypothesis parameter. Li is 
obtained by a Cholesky decomposition of the isotropic signal 

"^ \http: //astro, esa. int/Planck/science/performance/perf. top.htmI\ 
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covariance matrix d = (a^rnaj/^,). Using such an entropic 
prior instead of a Gaussian, as is the case when Wiener fil- 
tering, ensures a better reconstruction of non-Gaussian com- 
ponents such as the thermal and kinematic SZ effect. 

A maximisation of the posterior with respect to aim is 
therefore equivalent to minimising the function 



$(h^ 



X^i^im) - aS{Uim,m), 



(28) 



(29) 



since assuming Gaussian instrumental noise 

Pr{dim\h£m) OC CXp ["X^Cl^m)] , 

where x^ is the standard misfit statistic which utilises the 
isotropic noise covariance matrix N^ = (e^rnej,^,). 

To recover the SZ signal with the HSMEM, prior knowl- 
edge about cross-correlation between different components 
has not been assumed. Nevertheless, prior knowledge about 
the approximate shape and amplitude of the component 
power spectra is included in the analysis. This assumption 
is reasonable since preceding observations and data analysis 
of Planck observations themselves will provide prior con- 
straints on the component power spectra. 



4.2 Cluster detection algorithm 

To obtain the thermal SZ signal of single clusters the 
HSMEM recovered thermal SZ map is searched for local 
maxima above a certain threshold, which can be optimised 
to maximise the number of detected clusters and to minimise 
the number of false detections. A cluster is considered to be 
detected when the true position of its peak signal coincides 
with a local maximum found in the reconstructed map to 
within a defined off-set around the peak. In the work pre- 
sented here an off-set of one map pixel size (~ 2 arcmin; 
chosen to be just below the Nyquist sampling for the small- 
est Planck beam) has been allowed in the case the cluster's 
core radius angular extent is smaller than the size of a pixel. 
This applies to the vast majority of galaxy clusters. Other- 
wise the local maximum (in the following also referred to 
as SZ peak detection) has to lie within the core radius of 
the true cluster to be matched. This choice of off- set has 
been found to be optimal in order to obtain a high match- 
ing rate and avoid ambiguities arising from several maxima 
or clusters being matched to one of the other kind. Only 
maxima whose amplitude is greater than a fixed multiple 
of the estimated noise rms of the reconstructed map, anoise, 
(minimally Sxanoise) have been taken into account. The rms 
of the reconstructed map is obtained iteratively by masking 
maxima and their surroundings which are likely to belong 
to a cluster detection. Starting with the rms obtained over 
the whole map we exclude maxima above a threshold, for 
example 5 x anoise, and their environment which is assigned 
to them by flux integration. This process is repeated until 
the rms has converged. 

If a maximum can be matched to several clusters or 
several maxima can be matched to one cluster then an am- 
plitude and offset comparison is carried out in order to ob- 
tain the most likely match. In the case of extended clusters 
it can happen that due to noise (and background confu- 
sion) several maxima belong to the same cluster. Therefore, 
if maxima fall within each others estimated radius they are 
assumed to relate to the same extended cluster. 
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(d) (78 = 0.9 (xnorm) 



(e) (78 =0.7 (xnorm) 



(f) (78 = 1.0 (xnorm) 



Figure 2. 12.5 x 12.5deg^ thermal SZ realizations of the same patch of sky for three different values of cfq and two different M-T scaling 
relations (see equation Eland text) for the ACDM model (see TableQ used to make simulated Planck 5'wri'ey or observations. The value 
of (78 and the M-T relation used in the particular realization are given below each panel. The maps are shown in units of the frequency 
independent Compton y parameter. 



The integrated Comptonization parameter Y and clus- 
ter radius Vd associated with a SZ peak detection are esti- 
mated by spherical flux integration. At the radial extent at 
which the change in the interior integrated flux falls below 
a certain threshold, i.e. the contribution which is expected 
to be due to 3 x cFnoise^ the cluster is cut off. This ensures 
that the signal-to-noise ratio of cluster detections is at least 
3. 

Even though recently more sophisticated algorithms 
have been suggested, such as the pseudo optimal filter (Her- 
ranz et al. 2002) and MCMC methods (Hobson & McLach- 
lan 2003) , we decided to use the presented algorithm since it 
is efficient and has the ability to provide robust results (see 
section 15. 2f) . The advantage of the algorithm presented here 
is that, besides the assumption of the cluster to be spherical, 
the flux estimate is independent of the chosen cluster proflle. 
This is not the case for a pseudo optimal matched filter or 
MCMC application. The assumption that the cluster shape 
is spherical should be a fair approximation for most clusters. 



since X-ray observations show that the clusters' median el- 
lipticity is ^ 0.2 (see Mohr et al. 1995). 



We tested the assumptions on which the matching algo- 
rithm relies, such as the choice of the acceptance region and 
the association of multiple SZ peaks to a single extended 
cluster, and found, even though further optimisation may 
still be possible, that our choices are reasonable and do not 
affect the obtained completeness and purity estimates to a 
great extent. Nevertheless, problems which we found are dis- 
cussed in the text below (see section 15^ . Furthermore, we 
note that it is most likely that the final analysis performed 
on the Planck data will consist of a combination of several 
methods making optimal use of their different advantages, 
such as speed, detection efficiency and prior knowledge, at 
different levels of the analysis to maximise the number of 
cluster detections and the speed of the analysis. 
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Figure 3. The same as Figure |2l for the rCDM model (see Table CJ. Here a value of crs =0.6 has been used together with the different 
M-T relation normalisations to produce the thermal SZ simulations. 



5 RESULTS AND DISCUSSION 

Based on the modelling described above, full-sky maps of 
the SZ effects are constructed. Patches of the Compton y 
maps of the thermal SZ effect for different values of as and 
different M-T relations are shown in Figures |21 and El for the 
ACDM and rCDM cosmology respectively. The models have 
been chosen to cover a wide range of parameter space. In par- 
ticular, models which are supported by recent observations, 
such as the WMAP mission and X-ray cluster and weak lens- 
ing measurements, have been included. Moreover, the model 
choice has been influenced by observational data of the CMB 
at high angular multipoles (1500 < / < 3000). To visualise 
the differences occurring from changes in as and of the M-T 
relation, the same sky patch is shown for the ACDM and 
rCDM models respectively. Due to different structure for- 
mation histories the surface density of extended low redshift 
foreground clusters is larger in the rCDM model, whereas 
the ACDM models possess more high redshift background 
clusters. 

In Figures 31 and |5| the primordial CMB normalised to 
WMAP measurements at low multipoles and thermal SZ 
power spectra (in the following referred to as SZ power spec- 
tra) are plotted at 300 GHz. The included points with er- 
ror bars are observational results from ACBAR (Kuo et al. 
2002) and CBI (Mason et al. 2001). The measurements have 
been converted from the particular instrument observing fre- 
quency to 300 GHz after subtraction of the primordial CMB 
contribution by assuming that the remaining measured ex- 
cess is entirely due to the SZ effect and by using the spectral 
SZ frequency dependence (see equation|2|. The fraction con- 
tributed by the primordial CMB has been estimated from 
an extrapolation of the WMAP results best-fitting primor- 
dial CMB power spectrum to higher multipoles, which has 
been normalised to the height of the first acoustic peak. 
Further, it has been assumed that the primordial CMB esti- 
mate is free of error. Note that, at la level, no excess is re- 
quired to explain the ACBAR measurements. However, un- 
der the made assumptions the measured values of ACBAR 
and CBI are consistent with each other, the SZ interpret a- 
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Figure 4. The power spectra of the primordial CMB and the SZ 
effect for different values of ag and M-T relations (see panel and 
text) at 300 GHz for the ACDM simulations. While the upper 
panel shows a log-log plot, in the lower panel the SZ power spec- 
tra amplitudes are plotted on a linear scale in order to emphasise 
the different slopes and multipoles at which the SZ power spec- 
tra peak given different values of crs- The upper panel includes 
the CMB excess power spectrum amplitudes as measured by the 
instruments, CBI (circles) and ACBAR (squares). The measure- 
ments have been converted from their values at the particular 
instrument observing frequency to 300 GHz after subtraction of 
the primordial CMB contribution (see text for details) by assum- 
ing that they are entirely due to the SZ effect and by then using 
the spectral SZ frequency dependence (see equation |2|- 
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Figure 5. The same as Figure HI for the rCDM simulations. 



tion and the primordial CMB contribution estimated on the 
basis of the WMAP measurements. The consistency of the 
CBI and ACBAR excesses has also been found by Goldstein 
et al. (2003) by simultaneously fitting primordial CMB and 
SZ power spectra to the observations without using WMAP 
constraints. While the high-as ACDM model SZ power spec- 
tra fit the observed excess measurements well, the lower as 
ACDM models tend to lie below the excess measured by 
CBI and ACBAR. This finding is consistent with Bond et 
al. (2002) who obtained as = 1.0 using SZ power spectra 
from hydrodynamical simulations to fit the CBI excess for 
ACDM cosmologies. The rCDM as =0.6 SZ power spectra 
slightly overpredict the SZ power compared to the observed 
excess (but they are excluded by less than la). Neverthe- 
less, for the presented work it is important to include a wide 
range of models in order to study model dependent selection 
effects of the employed cluster detection algorithm. 



5.1 The SZ power spectrum relation to as 

Since the value of as has a significant effect on the abun- 
dance of galaxy clusters and thus on the measured power 
of SZ anisotropics, in the following the dependence of the 
thermal SZ power spectrum amplitude on as is investi- 
gated on the basis of the performed simulations. Predictions 
by various authors (Komatsu & Kitayama 1999; Seljak et 
al. 2001; Zhang & Pen 2001; Komatsu & Seljak 2002; Sadeh 
& Rephaeli 2004) expect Ci oc (Qbh) ag~ . Based on such 
predictions Bond et al. (2002) estimated as ~ 1.0 from the 
excess power at high / values measured by the CBI experi- 
ment using hydrodynamical simulations of the SZ effect. 

Most of the previous work on the SZ power spectrum 
dependence on as utilised either a standard Press- Schechter 
mass function (Press & Schechter 1976) or hydrodynamical 
simulations of small sky patches to obtain SZ power spectra. 
Komatsu & Seljak (2002) used the Jenkins mass function 



1 



Figure 6. The dependence of the power law index a on the multi- 
poles / for the virial (upper panel) and the M-T relation obtained 
from X-ray observations (lower panel) for our modelling of the SZ 
sky. 



(Jenkins et al. 2001), on which also the mass function of 
Evrard et al. (2002) is built, and a cluster pressure profile 
fitted to hydrodynamical simulations to evaluate the depen- 
dence of the SZ power spectrum coefficients on cosmological 
parameters, especially as, for the Poissonian contribution 
at small angular scales. Being based on the HV simulations 
our simulations ensure a realistic modelling of the clustering 
contribution of galaxy clusters to the SZ power spectrum. 

Moreover, the mass cut can affect the SZ power spec- 
trum estimation. For the simulated A and rCDM models we 
performed a mass cut at Md ^ 3 x 10^^ h~^ Mq . Komatsu & 
Seljak (2002) find that a mass cut as low as 5 x 10^^ h~^ Mq 
is sufficient to obtain the SZ power spectrum independently 
of mass boundaries for scales up to / ~ 10^. For such low 
'cluster' masses the applicability of the scaling relations as 
derived in sectional is questionable, since observations indi- 
cate a steepening of the M-T relation at mass scales of the 
order of 1 x 10^^/i~^ M© (Finoguenov et al. 2001; Borgani 
et al. 2002) . On the angular scales of interest the performed 
mass cut suffices to ensure independence of the SZ power 
spectrum. This has been tested by a comparison of the SZ 
power spectra obtained by performing different mass cuts 
and on this basis the performed mass cut has been found 
to be sufficient. For example, if we raise the mass cut from 
3 X W^h-^MQ to 5 X lO^^/i"^M0, a decrease of the SZ 
power spectrum amplitude at / = 5000 of approximately 10 
per cent is found. Anyway, since most of the clusters of the 
Planck survey will be unresolved (see section I53|) . Planck's 
ability to study substructures and cluster physics will be 
(very) limited. Therefore, small mass clusters are mainly in- 
cluded to provide a realistic SZ background noise. 

Nevertheless, we present in the following an estimation 
based on our cluster modelling of the dependence of the 
thermal SZ power spectrum amplitude on as over a wide 
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range of angular scales /. Expressing the relation between 
the SZ power spectrum amplitude and as by 



C^ocag", 



(30) 



a is obtained for two particular realisations of different as 

by 



«(/) = 



\og(Ci{a'i^)/Ci{af 
log (-^^V-f ^) 



(31) 



In Figure El the dependence of a on / is plotted for different 
combinations of cTg and cTg for the used M-T relations (see 
section 13.2. If) . 

Given our modelling of the SZ effect, a(l) is found to 
exhibit a significant negative slope. While for the ACDM 
models at low multipoles a adopts values of approximately 
6 to 7, on small angular scales (/ > 3000) it has fallen to 
values of at least a unit less. This behaviour can explain the 
wide range of values found in the literature for the power- law 
index a. Sadeh & Rephaeli (2004) obtained from analytical 
estimates of the SZ power spectrum using the cluster gas 
profile of equation (fTS|) Ci oc a^ at high /s, while previous 
studies at lower Is estimated Ci oc cTg"^ (see e.g. Komatsu 
& Kitayama 1999). Thus it has to be taken into account 
for which multipole range a is evaluated. Partially, this be- 
haviour is due to the clustering contribution which is most 
significant on large angular scales and has been shown to 
possess a stronger scaling with as (Komatsu & Kitayama 
1999). 

In the case of hydrodynamical simulations the contri- 
bution due to substructure, which is expected to increase 
the SZ power on small angular scales, possibly modifies the 
behaviour of the power-law index a at high multipoles. As 
some hydrodynamical simulations have been found to dis- 
agree with each other, it is at present impossible to make 
robust predictions. For example, the mesh-based GADGET 
code of Springel, White & Hernquist (2001) seems to predict 
more small-scale power than the simulations of da Silva et 
al. (2001) using HYDRA, a SPH code. Since the SPH sim- 
ulations of Bond et al. (2002) agree in shape with Springel, 
White & Hernquist (2001), the disagreement is unlikely to 
be dependent on the differences of the applied methods. The 
different shape of the high-/ SZ power spectra of various hy- 
drodynamical simulations could be due to a resolution ef- 
fect. The higher the resolution of the simulation the more 
substructure is revealed which increases the amplitude of 
the power spectrum at high multipoles and also flattens its 
peak, da Silva et al. (2004) state that they have not yet at- 
tempted a resolution study of their code. Anyway, since hy- 
drodynamical simulations are computationally expensive it 
is currently unfeasible to simulate volumes of the size needed 
for realistic full-sky simulations of Planck observations. 

Furthermore, in the case of the M-T relation derived 
from X-ray measurements the power-law index systemati- 
cally takes on smaller values than for the virial M-T rela- 
tion. This is due to the difference of the cluster Compton y 
profile induced by the change in temperature corresponding 
to a fixed cluster mass. Since the cluster radius and central 
electron gas density scale only with cluster mass and red- 
shift (see section El •, the y profile of the cluster is altered 
by the increase of Te . In our modelling the outskirts of the 
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Figure 7. Scatter plot of the reconstructed versus the actual 
cluster Comptonization parameter Y for clusters which have been 
matched to a SZ peak detection by the described algorithm (see 
section 221 for the fiducial ACDM model (see Table CJ. Based on 
the cluster detection algorithm the upper left, the upper right and 
the lower left panel show the scatter of the reconstructed cluster 
fluxes for clusters with a peak detection of 5, 4 and 3 times Onoise- 
The long dashed green line gives the optimal reconstruction. The 
horizontal red line represents a conservative survey Y limit esti- 
mate. It is chosen to lie well above the 'confusion plateau' level. 
The dotted vertical blue line shows the same flux limit for Yinput 
{Yinput ~ 2 X lO^'^arcmin^). Due to cluster-cluster confusion even 
low flux clusters of y ~ 1 x 10~^arcmin^ can have an associated 
Yrecon of scvcral 1 X 10~^arcmin^. The dashed red line gives a 
tolerant upper limit for when a cluster-SZ peak match is taken 
to be correct. This has been estimated from the skewness of the 
scatter around the optimal cluster flux reconstruction and is cho- 
sen to be 3 X Yinput (see text for detailed explanation). The lower 
right panel shows all detections. 

cluster are more affected by this modification of the profile 
than its core and therefore the effect is more significant at 
low multipoles. 

As Figure shows the power-law index a takes on a 
lower value in the case of the rCDM models. This agrees 
with the analytic expectations derived in Komatsu & Seljak 
(2002). Moreover, the peak of the SZ power spectrum is 
shifted towards smaller angular scales for smaller values of 
as (see lower panels of Figures 0] and |3 due to the rescaling 
of cluster masses (see equation ^^ and thus resizing of the 
clusters. 



5.2 The reconstructed SZ flux number count 

First we consider the concordance cosmological model, the 
ACDM model, which is supported by the WMAP measure- 
ments (Bennett et al. 2003), together with the virial M-T 
relation as our fiducial model. In particular the adopted 
cosmological parameters are Qm = 0.3, ^^a = 0.7, Hq = 
100/ikms"^Mpc"^ with h = 0.7, Qb = 0.04, n = 1 and 
as = 0.9 since these model parameters yield a reasonable 
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(xnorm) 


^cl (Yinput 
f^cl {Brecon 


> 5 X 
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10-3) 

10-3) 


3043 
3070 




4620 
4800 


ricl {Yinput 
^cl [j-recon 


> 4 X 
>4 X 


10-3) 
10-3) 


4134 
4194 




6243 
6602 


nd {Yinput 
^cl yj^recon 


> 3 X 
>3 X 


10-3) 
10-3) 


6030 
6200 




9264 
9803 


nd {Yinput 
^cl [j-recon 


>2 X 
> 2 X 


10-3) 
10-3) 


10441 
10523 




16051 
16729 


rici {Yinput 
^cl [^recon 


> 1 X 

> 1 X 


10-3) 
10-3) 


25410 
24160 




39689 
36155 



Table 4. The same as Tabled for the rCDM model. 



fit to the current strong cosmological constraints and thus 
provide a robust basis for making realistic predictions. Fur- 
thermore, these parameters agree with the ACDM HV N- 
body simulations used to obtain the thermal and kinematic 
SZ maps. 

Figure Qshows the scatter of the reconstructed SZ clus- 
ter integrated Comptonization parameters versus their real 
ones. The reconstructed ys have been obtained by apply- 
ing the algorithm as described in section 14.21 Only clus- 
ters whose real Comptonization parameter Y lies above 
1 X 10~^arcmin^ have been matched to detected SZ peaks. 
This assumed lower sensitivity limit of Y has been suggested 
by analytical considerations (see e.g. Bartelmann 2001). In 
Figure |7| the scatter of Y of detected clusters is shown for 
different peak detection thresholds (5, 4 and 3 x (Tnoise)- 
Even though the majority of detections with a peak above 4 
or ?>xanoise scatters symmetrically around the line of opti- 
mal detection and therefore can be believed to be matched 
up correctly, there is also a significant fraction which is 
scattered upwards to larger fluxes. These detections are ei- 
ther false matches or - which is more often the case - up- 
scattered due to cluster-cluster confusion. Thus for lower 
real cluster fluxes the reconstructed integrated Comptoniza- 
tion parameters level off. The shown survey limit estimate 
(Ysurveyiim — 2 X 10~^arcmin^) is chosen to lie well above 
this flux 'plateau' to ensure that the reconstructed flux dis- 
tribution represents closely the real one. As Figure Q illus- 
trates, the contribution of 3— and 4 — anoise peak detections 
above the sensitivity estimate for Planck is small. 

The reconstructed fluxes of some clusters are upscat- 
tered due to cluster- cluster confusion and even for the 
5 — (Tnoise dctcctious the scatter around the optimal re- 
construction is slightly skewed. Hence, there happen to be 
some outliers whose reconstructed Y is several times their 
Yinput- One way to deal with this is the introduction of a 
liberal upper y- limit that ensures that these questionable 
cluster reconstructions can be excluded from the sample. 
Thus by enforcement of this limit clusters with a recon- 
structed Yrecon abovc 3 X Yinput are regarded as false de- 
tections. On first sight this choice seems to be quite high, 
but one has to keep in mind that clusters confusing each 
other often have comparable masses and are situated next 
to each other in redshift space. In the following estimates 



of detection purity and completeness of the presented clus- 
ter detection method are given with and without taking 
account of this upper limit. Furthermore, these character- 
istics of the algorithm are shown for all cluster detections 
with associated peaks above 3 x (Tnoise as well as only for 
the 5 — (Tnoise pcak detections. Since the reconstructed in- 
tegrated Comptonization parameters of high-flux clusters 
scatter almost symmetrically around the line of optimal flux 
reconstruction, this indicates that the estimate of the total 
cluster flux obtained by our algorithm for clusters above 
Ysurveyiim — 2 X 10~'^arcmin^ is fairly unbiased. The level- 
ling off of the reconstructed ys at lower fluxes is mainly due 
to clusters being confused with each other. 

The completeness and purity obtained by the applied 
algorithm are plotted in Figure |H1 for the full-sky simula- 
tion and for a Galactic cut of Galactic latitude \b\ < 20°. 
The Galactic cut has been performed in order to evalu- 
ate the influence of Galactic residuals in the thermal SZ 
reconstruction on the quality of the obtained cluster num- 
ber count. For the full- sky cluster number count, the clus- 
ter survey is more than 83 per cent complete and the pu- 
rity exceeds 88 per cent above the estimated survey limit, 
Ysurveyiim = 2 X 10~'^arcmin^ , for all clusters without en- 
forcement of an upper matching limit. While the complete- 
ness is unaffected by the introduction of the proposed upper 
y- limit, the purity estimate is reduced by approximately 
5 per cent. This is the case when all detected clusters are 
considered as well as for only considering 5 — anoise peak 
detections. In any case neither the purity nor the complete- 
ness fall below 80 per cent for Ysurveyiim = 2 X 10~^arcmin^ 
for a cluster survey performed over the full sky. In the case 
of the Galactic cut, the completeness is increased by a few 
per cent, whereas the purity is barely improved for clus- 
ters with y above 2 x 10~'^arcmin^. Thus the Galaxy does 
not introduce significantly false detections into the recov- 
ered cluster sample. Some clusters are missed due to Galac- 
tic contamination. However, the improvement of the com- 
pleteness is mostly cancelled by the reduction in the num- 
ber of detected clusters which leads to an increase in the 
Poissonian sampling error of about the same percentage at 

^ surveylim ^ -^ X lU arcmiU . 

Note that the resulting completeness and purity also 
depend on the reliability of the cluster- SZ peak matching. 
Since the number of clusters and SZ peaks which have to 
be matched is large, the matching has to be automated. 
Even though an optimisation has been performed, it occa- 
sionally happens that the algorithm fails to match clusters 
and SZ peak detections correctly. The number of false de- 
tections and undetected clusters are both increased if a SZ 
peak which actually corresponds to a cluster is not matched 
up. This occurs, for example, in the case in which two clus- 
ters are confused with each other by beam-convolution and 
the SZ peak corresponding to both lies outside the match- 
ing acceptance region. Therefore, the quoted values have 
to be regarded as conservative lower limits of the cluster 
finding algorithm. This is supported by the comparison of 
the number of detected SZ peaks whose flux estimation ex- 
ceeds the sensitivity limit and the number of clusters above 
Ysurveyiim = 2 X 10~'^arcmin^ which are contained in the 
map (see Tables El and 2J . 

Figure EJshows the differential number redshift distribu- 
tion of clusters with and without a performed Galactic cut 
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Figure 8. In the left panel the same scatter as in FigureQis shown (for comparison), while a Galactic cut of ±20° Galactic latitude has 
been performed to obtain the scatter in the right panel. In this Figure the completeness and purity of the reconstructed cluster sample 
are shown. In the purity and completeness plots the red solid line represents the case in which all detections with a peak amplitude above 
3 X (Jnoise are taken into account, the blue long-dashed the case which utilises only detections whose peak amplitude is > 5 x dnoise- 
The black dotted line and the green short-dashed line show both cases when the matches which scatter above the red dashed diagonal 
line are not approved as detections. In the left panel the magenta short-dashed line gives the completeness of the survey when multiple 
clusters are allowed to be matched to a single detection peak. 
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0.7 (xnorm) 
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1.0 (xnorm) 
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10-3) 
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nci (Yinput > 4 X 

UcliYrecon > 4 X 


10-3) 
10-3) 
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379 
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1348/903 
1325/905 




2053 
2062 
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3186 


ncl (Yinput > 3 X 
T^clK^recon ^ <J X 


10-3) 
10-3) 


480 
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899 
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2083/1403 
2000/1371 
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3350 
3180 


4801 
4920 


rici {Yinput > 2 X 

ricliYrecon > 2 X 


10-3) 
10-3) 


930 

1478 


1741 
2360 


3781/2561 
3659/2532 




5770 
5998 


5941 

5784 


8701 
9035 


ncl (Yinput > 1 X 
ncl(Yrecon > 1 X 


10-3) 
10-3) 


2791 
5506 


5263 

7764 


10228/6735 
10744/7405 




16114 
16626 


15773 
15736 


23826 
23358 



Table 3. Cluster number counts above certain limits of the Comptonization parameter Y for the ACDM models with different ag and 
for the two utilised M-T relation normalisations. In each case the expected input and recovered number of clusters above the assumed 
Y limit (in arcmin^) is given. In the virial ag = 0.9 case the stated numbers are for a full-sky survey as well as for one with a Galactic 
cut of ±20° Galactic latitude (utilising 2/3 of the sky). 



as compared to the expected one for y > 2 x 10~^arcmin^. 
The largest difference between the expected distribution and 
the recovered ones occurs at low redshifts. This is mainly 
due to cluster- cluster confusion by beam-convolution. Due 
to Une-of-sight projection, an extended foreground cluster 
is more likely to overlap with background clusters. Part of 
this observed effect, however, is caused by the reconstruc- 
tion noise, since extended foreground clusters have smaller 
signal- to- noise peak detections than high-redshift clusters 



which appear more compact when projected along the line- 
of-sight. At higher redshifts the differential reconstructed 
number counts sometimes exceed those expected. This is 
due to clusters whose real fluxes are just below the survey 
threshold and being upscattered due to confusion of their 
SZ signal by projected overlapping clusters. The Galactic 
cut distribution yields slightly better results than the full- 
sky one. 

In Figure HTTl results of the application of the algorithm 
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Figure 9. The differential cluster number redshift distribution 
per square degree for the HV cluster catalogue (black solid line), 
for the full-sky recovered galaxy cluster survey (green dashed 
line) and when performing a Galactic cut of ±20° Galactic lati- 
tude (red dotted line) for an assumed survey sensitivity limit of 
^surveylim ^ 2 X 10 arcmm . 



to models with cosmological parameters which differ from 
the concordance model are shown. This has been done to 
test the algorithm which has been optimised for the fiducial 
ACDM cosmology. The models are chosen to incorporate 
a wide range of as which is still allowed by observations 
to ensure that in any case the method delivers reliable re- 
sults. The WMAP measurements on their own do not put 
strong constrains on as- Hence, as = 0.7 has been taken 
as the lower limit of the allowed range as this value is cur- 
rently supported by X-ray cluster survey analyses (Allen et 
al. 2003) . On the other hand weak lensing observations tend 
to favour higher values (see e.g. Hoekstra et al. 2002; Bacon 
et al. 2003). Thus as = 1.0 may be considered as an up- 
per limit of the allowed range. Besides the favoured ACDM 
model we include a flat model dominated by non-relativistic 
matter without a cosmological constant (rCDM model) in 
our investigations. 

We find, apart from the low-as ACDM model, the esti- 
mated completeness as well as purity obtained for the addi- 
tional models are of equivalent or even improved quality in 
comparison to estimates for the concordance model and the 
virial M-T relation. 

In the case of the low-as ACDM models noise intro- 
duced by the HSMEM itself and residual noise of other 
components become more significant due to the lower SZ 
signal amplitude. The completeness and purity for all clus- 
ter counts with Y above the survey limit are reduced to 56 
and 68 per cent respectively (without enforcement of an up- 
per matching limit). Thus, the applied algorithm does not 
quite suffice anymore to provide a reconstruction quality 
from which tight constraints on cosmological parameters. 



as and Qm, can be derived or which can establish a robust 
basis for successful follow-up observations. 

Moreover, for models with a larger variance of matter 
fluctuations on scales of 8 /i~^Mpc, the cluster confusion is 
more significant. Due to the mass rescaling obtained by em- 
ploying equation ((^ a cluster gets associated with a larger 
mass and therefore has a larger extent. This results in a more 
significant cluster overlap when projecting clusters SZ signal 
onto the two-dimensional full-sky map. As a result the com- 
pleteness slightly decreases due to increasing cluster-cluster 
confusion. On the other hand the purity is slightly increased 
due to the more dominant SZ signal. A similar effect is gen- 
erally found when the cluster temperatures are increased by 
a change of the M-T scaling relation normalisation. Then 
due to beam- convolution neighbouring clusters affect each 
other more strongly. Even though it can be argued that in 
any case of cluster-cluster confusion a confirmation proce- 
dure would be able to match the contributing clusters and 
even disentangle their contributions, these more optimistic 
considerations have not entered the presented results. 

Tables El and 2] give the detected cluster numbers above 
a declining threshold for each simulation. For comparison, 
the number of clusters whose flux exceeds the particular 
threshold and has been obtained directly from the cluster 
catalogues is also given. For Ysurveyiim = 2x 10~^arcmin^ we 
expect Planck to detect several 1000 clusters. If Ysurveyiim 
can be lowered to 1 x 10~^arcmin^, Planck will find a total 
number of clusters of the order lO'^. Since the recovered and 
input cluster number show good agreement even for Yum = 
1 X 10~^arcmin^, it is likely that, when the matching of SZ 
peak detections with clusters is even more optimised, a flux 
limit of 1 X 10~'^arcmin^ may be reached. 

So far we have taken into account the quality of the clus- 
ter flux recovery to estimate a sensitivity limit for a Planck 
cluster survey. Previous work often considered a cluster to 
be detected when it could be successfully matched with a 
peak in the reconstructed SZ map and allowed even several 
clusters to be matched with a single SZ peak. It can be ar- 
gued that follow-up observations at higher resolution or at 
other wavelengths are able to disentangle the contributions 
of several clusters to an extended SZ feature and to confirm 
Planck cluster detections. Such follow-up observations can 
be performed, for example, by ground-based SZ telescopes, 
such as the interferometers AMI (Kneissl et al. 2001), Amiba 
(Lo et el. 2000) and SZA (Mohr et al. 2002) in the north- 
ern hemisphere and the upcoming instruments in the south, 
such as SPT^ and ACT^. If we assume that the contributions 
of individual clusters can be separated, the completeness is 
significantly increased (see Figure Q). The obtained cluster 
sample is 83 per cent complete at Y^^rn = 1 x lO^'^arcmin^ 
and still has a purity of 78 per cent. 

For the fiducial ACDM model. Figure ITTl compares the 
effect of these assumptions on the cluster sample obtained 
for the two considered sensitivity limits by the differential 
redshift number distributions. It is found that due to cluster- 
cluster confusion, especially at higher redshifts {z > 0.3), at 
which the angular separation of clusters situated next to 
each other in redshift space is reduced, the recovered dis- 



^ http://astro.uchicago.edu/spt 

^ http://www.hep.upenn.edu/~angelica/act/act.html 
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(a) ACDM: erg = 0.7 (virial) 



(b) ACDM: erg = 1.0 (virial) 



(c) rCDM: erg = 0.6 (virial) 
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Figure 10. The scatter, completeness and purity for all performed cosmological SZ realisations apart from the fiducial ACDM one are 
shown. 



tributions using the recovered cluster fluxes exceed the ex- 
pected cluster numbers significantly for both survey limits. 
This is caused by an increased likelihood of clusters with 
fluxes below the particular limit to be included in the sam- 
ple, since all clusters contributing to the connected peaked 



SZ feature are taken into account. On the basis of Planck 
observations analysed with the presented cluster detection 
method, it is not possible to separate clusters contributing 
to a single SZ peak detection. Thus, follow-up observations 
are unavoidable in order to disentangle clusters and remove 
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Figure 11. The differential cluster number redshift distribution 
per square degree for the HV input cluster catalogue and for the 
full-sky galaxy cluster detections under the assumption that a 
cluster is detected when its matching region includes a recon- 
structed SZ peak. Thus, multiple clusters can be matched to a 
single peak. In the case of the peak-detected clusters above a given 
flux survey limit the distributions are shown for the true cluster 
flux as well as for the flux associated with the peak matched to the 
cluster (different distributions are indicated by different colour 
and line-styles). The plot shows the distributions for assumed sur- 
vey sensitivity limits of 2 x 10 ~ 



low-flux clusters contributing to a SZ peak detection above 
the flux limit from the sample. Then it is feasible to reach 
the distributions which are given for the cluster peak de- 
tections assuming knowledge of the true cluster fluxes (see 
Figure ITT)) . Further, Figure ITTl shows that the recovered dis- 
tributions (using cluster input as well as recovered fluxes) of 
the higher survey flux limit (y = 2 x 10~^arcmin^) resem- 
ble the true input distribution more closely than it is the 
case for the lower flux limit, Yum = 1 x 10~^arcmin^. In the 
following the more conservative estimate of the sensitivity 
limit, Ysurveyiim = 2 X 10~'^arcmin^ , is assumed and only 
single cluster-peak matches are allowed. 

Based on the reconstructed number counts, the purity 
and completeness estimate, we can examine how likely mod- 
els exclude each other given the count of a certain model and 
comparing it to the expected count of the ones that yield 
similar cluster numbers above the conservative survey limit. 
This is flrst done by assuming that our estimates of the per- 
centage of impure and missing clusters are standard devia- 
tions of a Gaussian distribution centred on the reconstructed 
and expected cluster number respectively. Furthermore, we 
account for the cosmic variance in the number of clusters. 
Assuming a M-T relation, we flnd that given the recovered 
count of a particular model, each of the other models consid- 
ered here can be excluded by at least 1.7a. This minimal ex- 
clusion is obtained in the case of the reconstructed ag = 0.9 
ACDM count when compared with the expected count of 



the as = 1.0 ACDM model. However, the ag = 0.6 rCDM 
cluster count excludes the as = 1.0 ACDM model by 2.6a. 
When we take our estimates of impurity and incomplete- 
ness to be absolute values, which have an error themselves, 
then the probability of exclusion of the respective models 
with similar expected number counts exceed 5a in all cases. 
Since these estimates have been obtained by a large number 
of single cluster detections and we can regard them as the 
probability that a cluster above the survey limit is detected 
or a real detection respectively, this assumption should be 
valid and we can assume a Poissonian sampling error in the 
completeness and purity. Nevertheless, we first considered 
the more conservative case of the estimated incompleteness 
and impurity to be the first moments of a Gaussian distri- 
bution of the expected and recovered number count. 

As expected the cluster number exceeding a given fiux 
limit is increased for the M-T scaling relation normalised by 
X-ray observations (xnorm) . An uncertainty in the normali- 
sation of the M-T relation introduces extra degeneracies be- 
sides the ones present between cosmological parameters. For 
example, the cluster number yield of the ACDM ag — 0.9 
model employing the xnorm M-T scaling relation is compa- 
rable to the one obtained for the ACDM as = 1.0 cosmology 
using the virial M-T relation. 

5.3 Cluster size distribution and cluster resolution 

To study galaxy cluster physics resolution is important. For 
the Planck satellite, the large beam sizes mean that most 
clusters will be unresolved. Nevertheless, a minor fraction of 
the cluster sample, nearby massive clusters, will be resolved. 
Besides giving an estimate of the expected percentage of re- 
solved clusters contained in the sample above the survey 
limit estimate {Ysurveyiim — 2 X 10~^arcmin^), we further 
investigate how well cluster radii are recovered by the ap- 
plied cluster finding algorithm. 

In Figure IT^ we plot the distribution of the cluster an- 
gular sizes for a plurality of cluster size definitions. First 
of all the distribution of the virial radii obtained from the 
scaling relation (|21f) is shown. Furthermore, the radius dis- 
tribution of the beam- convolved cluster extent at which the 
cluster SZ signal has dropped to half of its peak amplitude is 
given. This can be compared with the FWHM of the beam 
estimate of 9.7 arcmin for the HSMEM reconstructed SZ 
map. This beam estimate of the reconstructed thermal SZ 
map has been obtained by applying a x^-^i of the convolved 
input map to the reconstructed thermal SZ map. If we de- 
fine that a cluster is convolved in the case that this FWHM 
radius estimate exceeds the FWHM of the reconstruction 
beam, which ensures that a cluster significantly contributes 
to several beams, a fraction of approximately 10 per cent of 
the recovered clusters of fiux y > 2 x 10~^arcmin^ is found 
to be resolved. 

Moreover, for comparison the expected and recovered 
cluster radius distributions for the cluster extent at which 
the change of the integrated fiux with radius becomes equiv- 
alent to the contribution expected from 3 x anoise are shown 
in Figure IT^ Both distributions show the same onset at low 
cluster radii (approximately 8 to 9 arcmin). Due to over- 
estimation of the cluster radii by cluster-cluster confusion 
and noise, the distribution obtained from the reconstructed 
SZ map appears to be broader and less peaked. However, 
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Figure 12. The angular size distribution of clusters in the fiducial 
ACDM survey with Y >2x 10~^arcmin^. The solid black curve 
represents the distribution of the virial angular radius, the orange 
dotted-dashed line the distribution of the radius at half the am- 
plitude of the cluster SZ signal after a convolution with a beam of 
FWHM 9.7 arcmin. The distributions of the reconstructed cluster 
radii and the expected radius reconstructions neglecting cluster- 
cluster confusion are given by the green dashed and red dotted 
line respectively. The cluster is cut off at the radius at which the 
change of the integrated fiux with radius falls below the contri- 
bution expected to be due to 3 x (Jnoise- The blue long-dashed 
vertical line represents the FWHM angular extent of the beam 
estimate for the reconstructed map. 



for large extended clusters the reconstructed distribution 
again approaches the estimated one and the algorithm pro- 
vides thus a reasonable estimation of the cluster radii. This 
occasional overestimation of projected small and mid-sized 
cluster radii partially accounts for the upscattering of the 
recovered flux distribution. 



5.4 The reconstructed SZ poAver spectrum 

In addition to the cluster number count, the SZ power spec- 
trum also strongly depends on cosmology and can thus be 
used to constrain cosmological parameters. Therefore, it is 
worthwhile investigating how well the SZ power spectra of 
the different models can be reconstructed by the HSMEM. 
In Figures QSl and El we plot the HSMEM reconstructed 
SZ power spectra for each of the considered models. For 
comparison also the input as well as the beam- convolved 
SZ power spectra are shown. To obtain the beam-smoothed 
power spectra, a beam of FWHM of 9.7 arcmin has been 
assumed, since it provides the best x^-^i to the recon- 
struction of our fiducial ACDM model. Even though the 
HSMEM reconstructed power spectrum is intrinsically bi- 
ased (see Stolyarov et al. 2002), in the case of the fiducial 
ACDM model, the reconstructed SZ power spectrum agrees 
well with the beam- convolved input SZ power spectrum over 
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Figure 13. Comparison of the reconstructed thermal SZ effect 
power spectra in the case of the ACDM model {Qm = 0.3, A = 
0.7, h = 0.7; Table nj for different values of ag. The values of 
ag, which refer to the particular input and reconstructed power 
spectra, are indicated by the colour and line style. In addition 
we plot the input (thick lines) and beam convolved (thin lines) 
SZ power spectra. To obtain the beam convolved ones we assume 
a Gaussian beam of FWHM of 9.7 arcmin. For comparison the 
reconstructed primordial CMB power spectrum is also shown. 
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Figure 14. The same asHalfor the rCDM model {Qm = 1, A = 0, 
h = 0.5; TableCJ. 
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a wide range of multipoles. Only on medium angular scales 
(500 < / < 1000) the relative difference can be as large as 50 
per cent. Nevertheless, at lower multipoles the reconstructed 
power spectrum is in good agreement with the correspond- 
ing true one for most of the models. Only in the case of the 
low-cTs ACDM model (especially when using the virial M- 
T relation) the reconstructed SZ power spectrum appears 
to be biased. In particular, the reconstruction overpredicts 
the SZ power on large scales. This is due to reconstruction 
and residual noise of other components present in the recon- 
structed SZ map. Since the SZ power of the low-as ACDM 
models is approximately an order of magnitude smaller than 
for the high-cTs ACDM models, the obtained power spectra 
of the former models are more affected by reconstruction 
and residual noise. Hence, one obtains for these models an 
overestimated SZ power amplitude. 

Furthermore, as the comparison of the beam- convolved 
with the reconstructed power spectrum shows, models with 
larger amplitudes than the fiducial ACDM model repro- 
duce the true SZ power spectrum better at higher multi- 
poles and seem to be less convolved. This effect arises due 
to the relative increase of the SZ signal in comparison to the 
instrumental noise and other components with fixed am- 
plitude, such as the Galactic dust, in the high-resolution 
high-frequency channels of Planck. The ones with lower am- 
plitudes show the opposite behaviour. However, up to mul- 
tipoles of several hundreds the reconstructed power spec- 
tra are not significantly affected by the effective convolving 
beam of the reconstructed map. Thus given adequate prior 
knowledge about the spatial correlations of galaxy clusters, 
the outer cluster profile and the M-T relation, our simula- 
tions suggest that the reconstructed SZ power spectrum at 
low multipoles can provide constraints on cosmological pa- 
rameters in addition to the cluster number count. As shown 
by Komatsu & Seljak (2002) the constraints on the cos- 
mological parameters obtained by the SZ power spectrum 
can be used to break degeneracies between parameter con- 
straints deduced from cluster number counts. 



6 SUMMARY AND CONCLUSIONS 

The Planck Surveyor satellite offers the opportunity to com- 
pile a large sample of galaxy clusters over the full sky. This 
additional aim of the Planck mission - besides determin- 
ing the primordial CMB power spectrum at high angular 
resolution - can be a very powerful cosmological tool. We 
have performed full-sky multi-component simulations of fu- 
ture Planck SZ observations with a high degree of realism. 
This is the first time that on the basis of full-sky simulations 
and an application of a particular detection algorithm the 
SZ cluster gain of the Planck mission has been investigated. 
Based on cluster mass and redshift catalogues obtained 
from N-body simulations and a mass function we derived a 
model sky, using further cluster properties from scaling rela- 
tions. Our modelling utilised two different normalisations of 
the M-T scaling relation to explore the dependence of our re- 
sults on its choice. The one derived from X-ray observations 
is more likely to include cluster physics, such as preheat- 
ing and cooling, and gives a higher gas temperature for a 
given cluster mass. The theoretically derived one has been 
normalised with adiabatic hydrodynamical simulations and 



presents a lower bound. To model the extent of the cluster 
electron gas we used the /3-profile which is sufficient given 
the resolution of Planck. On the basis of this modelling we 
investigated the dependence of the SZ power spectra on as. 
We found that the power law index a{l) of the Ci depen- 
dence on as possesses a negative slope in /-space, which we 
argue can explain the wide range of estimates for a obtained 
previously by various authors. Moreover, all confusing extra- 
galactic and Galactic components known to be significant 
have been incorporated in our full-sky simulations of the 
Planck channel detections. By the current state of art our 
modelling is sufficient to yield reliable results. 

Further, we have applied a SZ cluster detection algo- 
rithm existing of two steps, the HSMEM to separate the 
SZ signal from other contaminating components and a clus- 
ter finder. In the present work the simulated performance 
of Planck as a SZ cluster survey instrument in combina- 
tion with our detection algorithm has been tested for sev- 
eral models including the ones favoured by current data. 
On the basis of the fiducial ACDM cosmology (see Ta- 
ble we estimated a conservative flux survey limit of 
Y — 2x 10~^arcmin^. Given this limit a sufficient complete- 
ness and purity is obtained for the majority of the investi- 
gated models. Only in the case of the as = 0.7 ACDM model 
the completeness and purity have been found to be signif- 
icantly lower than 80% for Ysurveyiim — 2 x 10~^arcmin^. 
In the case that cluster contributions to a single SZ peak 
detection can be disentangled, it has been shown that the 
survey limit can be lowered to 1 x lO^'^arcmin^. Our inves- 
tigations indicate that the cluster selection of the detection 
algorithm is non-trivial and shows a cosmological depen- 
dence mainly due to the strong cosmological variation of 
the SZ signal magnitude. Ideally, the cluster detection al- 
gorithm returns a cluster sample whose completeness and 
purity are cosmology independent, since otherwise one has 
to rely on prior information and mock simulations as per- 
formed in the present work. However, in practice this ideal 
case can hardly be reached. Nevertheless, prior constraints 
can be obtained from the power spectrum of the SZ effect 
and thus allow a self- consistent cosmological parameter es- 
timation based on the Planck SZ detection. Therefore, as 
suggested by our simulations besides serving as a basis for 
follow-up requiring a major observational effort the recov- 
ered Planck cluster catalogue itself can already be used to 
put constraints on cosmological parameters. A simple con- 
servative investigation based on our SZ detection algorithm 
in section 15.21 showed that a Planck SZ cluster survey will 
be able to distinguish between the considered vanilla models 
by at least a few a. Given our algorithm we will investigate 
how tightly cosmological parameters can be constrained by 
a Planck SZ survey in a future paper. 

Thus, as our results suggest Planck will not only be a 
powerful tool to constrain cosmology via primordial CMB 
fluctuations but will also have the ability to obtain indepen- 
dently constraints from the secondary effect of SZ distor- 
tions. 



7 ACKNOWLEDGEMENTS 

JG acknowledges support by an Isaac Newton Fellowship of 
the Isaac Newton Trust of the University of Cambridge. We 



© 0000 RAS, MNRAS 000, 000-000 



Sunyaev-ZeVdovich Cluster Survey Simulations for Planck 19 



especially thank Vlad Stolyarov for providing the HSMEM 
code and for several helpful discussions about its appli- 
cation to full-sky simulations. Moreover, we acknowledge 
the HEALPix collaboration for providing the pixelisation 
scheme used in this work and the Virgo Consortium for 
making the Hubble Volume N-body simulations publically 
available. 



REFERENCES 

Allen S. W., Schmidt R. W., Fabian A. C, 2001, MNRAS, 328, 
37 

Allen S. W., Schmidt R. W., Fabian A. C. &; Ebehng H., 2003, 
MNRAS, 342, 287 

Arnaud M., Pratt G. W., Pointecouteau E., Memorie della Soci- 
eta' Astronomica Italiana, the Proceedings of the EPIC Con- 
sortium (held on Oct 14-16, 2003 in Palermo) 

Bacon D.J., Massey R.J., Refregier A.R., Ellis R.S., 2003, MN- 
RAS, 344, 673 

Bartelmann M., 2001, A&A, 370, 754 

Bennett C.L., et al., 2003, ApJS, 148, 1 

Bennett C.L., et al., 2003, ApJS, 148, 97 

Birkinshaw M., 1999, Phys. Rep., 310, 98 

Plain A. W., Ivison R. J. & Smail I., 1998, MNRAS, 296L, 29B 

Bond J. R., et al., 2002, ApJ, submitted lastro- ph/0205386 I 

Borgani S., Governato F., Wadsley J., Menci N., Tozzi P., Quinn 
T., Stadel J., Lake C, 2002, MNRAS, 336, 409 

Bouchet F.R., Gispert R., 1999, NewA, 4, 443 

Burigana C, Danese L., De Zotti G., Franceschini A., Mazzei P., 
Toffolatti L., 1997, MNRAS, 287, L17 

Carlstrom J.E., Holder G.P, Reese E.D., 2002, ARAA, 40, 643 

Chalhnor A., Lasenby A., 1998, ApJ, 499, 1 

Danese L., De Zotti G., Francheschini A., Toffolatti L., 1987, ApJ, 
318, L15 

da Silva, A. C, Kay S. T., Liddle, A. R., Thomas, P. A., Pearce, 
F. R. & Barbosa, D., 2001, ApJL, 561, 15 

da Silva A.C., Kay S.T., Liddle A.R., Thomas PA., 2004, MN- 
RAS, 348, 1401 

Diego J.M., Vielva P., Martinez-Gonzalez E., Silk J., Sanz J.L., 
2002, MNRAS, 336, 1351 

Eke V. R., Cole S., Frenk C. S., 1996, MNRAS, 282, 263 

Eke V. R., Navarro J. F., Frenk C. S., 1998, ApJ, 503, 569 

Ensshn T. A., Hansen S. H., 2004 (|astro-ph/0401373 ) 

Ettori S., Tozzi P., Borgani S., Rosati P., 2004, A&A, 417, 13 

Evrard A.E., et al., 2002, ApJ, 573, 7 

Finoguenov A., Reiprich T.H., Bohringer H., 2001, A&A, 368, 
749 

Francheschini A., Mazzei P., De Zotti G., Danese L., 1994, ApJ, 
427, 140 

Frenk C.S., et al., 2000 ( astro-ph/00073^ 

Goldstein ,et al., 2003, ApJ, 599, 773 

GuU S. F., Skilhng J., 1990, The MEMSYS5 User's Manual. Max- 
imum Entropy Data Consultants Ltd, Royston 

Haslam C.G.T., Salter C.J., Stoffel H., Wilson W.E., 1982, A&AS, 
47, 1 

Herranz D, Sanz J.L., Barreiro R.B., Martinez-Gonzalez E., 2002, 
ApJ, 580, 610 

Hobson M.P., Barreiro. R.B., Toffolatti, L., Lasenby A.N., Sanz 
J.L., Jones A.W., Bouchet F.R., 1999, MNRAS. 306, 232 

Hobson M.P, Jones A.W., Lasenby A.N., Bouchet F.R., 1998, 
MNRAS, 300, 1 

Hobson M.P, Lasenby A.N., 1998, MNRAS, 298, 905 

Hobson M.P, McLachlan C, 2003, MNRAS, 338, 765 

Hoekstra H., Yee H.K.C., Gladders, M.D., 2002, ApJ, 577, 595 

Holland W.S., et al. 1998, Nature, 392, 788 

Hu W., Kravtsov A.V., 2003, ApJ, 584, 702 



Jenkins A., Frenk C.S., White S.D.M., Colberg J.M., Cole S., 

Evrard A.E., Yoshida N., 2001, MNRAS, 321, 371 
Jonas J.L., Baart E.E., Nicolson G.D., 1998, MNRAS, 297, 977 
Kaiser, N., 1986, MNRAS, 222, 323 

Kay S.T., Liddle A.R., Thomas P., 2001, MNRAS, 325, 835 
Kneissl R., Jones M.E., Saunders R., Eke V.R., Lasenby A.N., 

Grainge K., Cotter G., 2001, MNRAS, 328, 783 
Komatsu, E. &; Kitayama, T., 1999, ApJL, 526, LI 
Komatsu, E. & Seljak, U., 2002, MNRAS, 336, 1256 
Kuo C.-L., et al., 2003, AAS, 34, 1324 
Lacey C, Cole S., 1993, MNRAS, 262, 627 
Lewis A., Chalhnor A., Lasenby A., 2000, ApJ, 538, 473 
Lo K.Y., Chiueh T., Martin R.N., 2000, AAS, 197, 5603 
Mason B., CBI Collaboration, 2001, AAS, 199, 3401 
Mather J. C, Fixsen D. J., Shafer R. A., Mosier C, Wilkinson 

D. T., 1999, ApJ, 512, 511 
Mohr J. J. et al., 1995, ApJ, 447, 8 

Mohr J.J., Carlstrom J.E., The Sza Collaboration, 2002, AMiBA 
2001: High-Z Clusters, Missing Baryons, and CMB Polariza- 
tion, ASP Conference Proceedings, Vol. 257. Edited by Lin- 
Wen Chen, Chung-Pei Ma, Kin- Wang Ng, and Ue-Li Pen. 
Navarro J. F., Frenk C. S., White S. D. M., 1997, ApJ, 490, 493 
Pointecouteau E., Arnaud M., Kaastra J., de Plaa J., 2004, A&A, 

accepted 
Pointecouteau E., Giard M., Barret D., 1998, A&A, 336, 44 
Pierpaoh E., Scott D., White M., 2001, MNRAS, 325, 77 
Press W. H., Schechter P., 1974, ApJ, 187, 425 
Reich P., Reich W., 1986, A&AS, 63, 205 
Rephaeh Y., 1995, ARA&A, 33, 541 
Sadeh S., Rephaeh Y., 2004, NewA, 9, 373 
Sarazin C. L., 1986, Rev. Mod. Phys., 58, 1 
Schlegel D.J., Finkbeiner D.P, Davis M., 1998, ApJ, 500, 525 
Sheth R.K., Tormen G., 1999, MNRAS, 308, 119 
Seljak U., BurweU J., Pen U.-L., 2001, PhRvD, 63, 3001 
Smail I., Ivison R. J., Plain A.W., Kneib J.-P, 1998, AAS, 192, 

4813 
Springel, V., White, M. & Hernquist, L., 2001, ApJ, 549, 681 
Stolyarov V., Hobson M.P., Ashdown M.A.J. , Lasenby A.N., 

2002, MNRAS, 336, 97 
Stolyarov V., Hobson M.P., Lase nby A .N., Barreiro R.B., 2004, 

MNRAS, submitted rastFo-ph70405494 ) 
Sunyaev R.A., Zel'dovich Ya. B., 1972, Comm. Astrophys. Space 

Phys., 4, 173 
Sunyaev R.A., Zel'dovich Ya. B., 1980, ARA&A, 18, 537 
Toffolatti L., Argiieso Gomez F., De Zotti G., Mazzei P., Franch- 
eschini A., Danese L., Burigana C, 1998, MNRAS., 297, 117 
White, M. 2003, ApJ, 597, 650 
Zhang, P.J. & Pen, U.L., 2001, ApJ, 549, 18 



© 0000 RAS, MNRAS 000, 000-000 



